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Control of intramolecular spin alignment is studied theoretically in a model of polyene-based 
molecular magnets in which delocalized n electrons are coupled with localized radical spins. In a 
previous paper [Phys. Rev. Lett. 90 (2003) 207203], we have demonstrated that charge doping is an 
effective way to realize controllable high-spin in the 7r-conjugated molecular magnets. In this paper, 
we clarify the dependence of spin-alignment on the exchange interaction between the localized spin 
and 7r electron and the electron-electron interactions. The antiferromagnetic exchange interaction 
' plays a role different from the ferromagnetic counterpart in doped molecules. To understand complex 

^ \ interplay of charge and spin degrees of freedom in the doped systems, we carry out a systematic 

. study on the phase diagram of spin alignment in the parameter space. The mechanism of the 

spin alignment is discussed based on the spin densities of n electrons. The calculated results are 
consistent with experiments, providing a theoretical basis for the control of spin alignment. 
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I. INTRODUCTION 

Since the discovery of high-spin ground state in 
dicarbene^^ molecular magnetism has been exten- 
sively investigated both experimentally and theoret- 
ically for decades A new class of purely or- 
ganic molecular magnets has received increasing at- 
tention as their intramolecular spin alignment is con- 
trollable by external stimuli such as charge doping 
or photoexcitation) 7 i 8 i 9 i 10 i 1:L i 12 i 13 i 14 i 15 i 16 as schematically 
shown in Fig. ^ Those molecules share a common fea- 
ture in structure: two stable radicals carrying unpaired 
electrons coupled with a 7r-conjugated moiety through 
exchange interaction. The two radicals are barely cou- 
pled to each other directly, but interact indirectly via ir 
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FIG. 1: Schematic picture of control of intramolecular spin 
alignment by charge doping or photoexcitation. 



can be controlled by charge doping in thianthrene 
bis(nitronyl nitroxide)^ The molecule is spin singlet 
(S = 0) in neutral state with very weak antiferromag- 
netic exchange coupling (—1.4 K) between the nitronyl 
nitroxide groups. Upon one-electron oxidation, i.e. one- 
hole doping, spin quartet (S = 3/2) ESR signals are 
detected in the temperature range of 5 — 150 K, in- 
dicating a low-spin to high-spin transition induced by 
doping. Optically controllable spin alignment has been 
demonstrated by Teki et al. in a molecule composed of 
dangling iminonitroxide radicals and 7r-conjugated moi- 
ety of diphenylanthraceneiiSiii The ground state of the 
molecule is determined to be a low-spin state, corre- 
sponding to weak antiferromagnetic coupling (—5.8 K) 
between the radical spins. After the molecule is ex- 
cited in low temperature (20 — 40 K) by short pulses 
of laser, spin quintet (S = 2) species are observed in 
10 /its by using the time-resolved ESR spectroscopy. This 
metastable spin quintet results from photoinduced par- 
allel spin alignment among the two stable radicals and a 
photogenerated triplet ir radical. Since these molecular 
magnets are purely organic, there is a great advantage 
of various chemical modifications onto the 7r-conjugated 
moiety. It might allow such molecular magnets to be 
building blocks of photoinduced bulk magnets similar to 
Prussian-blue analogs i^Z^ 

Many theoretical works have been carried out to 
propose and understand the mechanism of spin align- 
ment in molecular magneta 19 i 2 VV 2 i 23 i 24 i 2 'V 6 i 2 7i 28 i 2 V n . 
By a weak electron-correlation approach, Mataga pre- 
dicted that ferromagnetic spin alignment appears in sev- 
eral hypothetical hydrocarbons with 7r-conjugated elec- 
tron systems due to the topological nature of molecu- 
lar orbitals. 19 Ovchinnikov accounted for the net spins 
of planer alternate hydrocarbons by a strong electron- 
correlation approach^ These two approaches were later 
generalized to rigorous theorems in the Hubbard model 



2 



by Lieb^i Thus the dominant mechanism of intramolec- 
ular spin alignment in the ground state of 7r-conjugated 
molecules has been well established by those works: The 
alignment is governed by the topological rule based on 
the dynamical antiferromagnetic spin polarization effect 
of 7r electrons with on-site Coulomb repulsion. 

The topological rule can be applied only to the half- 
filled ground states, but spin alignment in doped or ex- 
cited molecules remains unsolved for general cases. In 
the strong-correlation limit with an infinite value of on- 
site Coulomb repulsion in the Hubbard model, Nagaoka 
proved rigorously that the ferromagnetic ordering is fa- 
vorable in the case of single-electron addition to or re- 
moval from half-fillingi^ Most discussions on the spin- 
alignment control so far have been given on the basis of 
individual molecular orbitals&ii* 1 ^ For high-spin states 
induced by charge doping, Sakurai et al. interpreted par- 
allel alignment of two radical spins in terms of the elec- 
tronic structures calculated by the unrestricted Hartree- 
Fock theory^ A mechanism of photoexcited high-spin 
states was discussed by Teki et al. using ab initio molec- 
ular orbital calculations based on the density functional 
theory 11 ! 13 

Approaches by model Hamiltonians have been widely 
used to study the neutral ground state of quasi-one- 
dimensional (Id) organic ferromagnets. Nasu2£ studied 
spin alignment in an idealized poly(m-diphenylcarbene) 
using a periodic Kondo-Hubbard model. Fang et alJtL 
and Chen et al~ studied a periodic Kondo Su-Schrieffer- 
Hegger(SSH) model 33 taking into account the strong 
electron-phonon interaction as well as the electron- 
electron correlation. In the previous papers^*^, we in- 
troduced a Kondo-Pcicrls-Hubbard model to study spin 
alignment of 7r-conjugated molecular magnets. It has 
been demonstrated that a spin singlet to quartet tran- 
sition can be induced by electronic doping into the it 
moiety, which is ferromagnetically coupled with local- 
ized spins. The intramolecular spin alignment depends 
sensitively on the topological structure of the molecular 
system. The doped high-spin state can be realized in ap- 
propriate molecular structures. The topological effect in 
the doped case is very different from that in the half-filled 
case. 

The aim of the present paper is to clarify the depen- 
dence of spin-alignment on spin exchange and electron- 
electron interactions, thereby giving an integral picture of 
spin-alignment control in polyene-bascd molecular mag- 
nets. The paper is organized as follows: In Sec. [Q] we in- 
troduce the model and the calculation method. The cal- 
culated results are discussed in Sec. IIIII We first study 
the spin alignment and the lattice deformation pattern 
with antiferromagnetic coupling as well as the ferromag- 
netic one in Sec. IIII Al The dependence of spin alignment 
on the molecular structure is discussed in Sec. IIIIBI To 
study the complex interaction between charge and spin 
in doped quantum systems, we discuss the dependence of 
spin alignment on parameters J and U in Sec. IIII CI All 
the results are summarized in Sec. IIVI 



II. KONDO-PEIERLS-HUBBARD MODEL 

In this paper, we use the same Kondo-Peierls-Hubbard 
model as in our previous papers^ 4,35 : a Peierls-Hubbard 
model on a A-site linear chain interacting with two local- 
ized quantum spins in the form of the Kondo coupling. 
The open-boundary linear chain corresponds to polyene, 
and its n electrons are subject to electron-lattice inter- 
actions in the form of the SSH coupling^ The localized 
spins correspond to the unpaired electrons of stable rad- 
ical groups attached to the main chain, as realized in 
some substituted polyacetylenesi 36 ! 37 The polyene-based 
molecular magnets described by this model is currently 
hypothetical but there is a possibility to synthesize them 
in future. 

Thus, our model is presented in the Hamiltonian, 
H = -J2ti(Ct s C i+1 , s +h.c.) 

+ U n ht n i,l + _ 91+1 ^ 

i i 

-J(S n ■ S T i + S i2 • S T2 ) , (1) 

U = t a + a(q t - qi+t), (2) 

Si = \c}aC t . (3) 

In Eq. (1), C\ s {Ci tS ) creates (annihilates) an electron 
with spin s(=| or J,) on site i. The transfer integral ti 
is introduced for the nearest-neighbor sites, and is com- 
posed of a constant magnitude to, and a variable one 
that depends on the displacement qi of site i with the 
SSH coupling constant a. On-site Coulomb repulsion 
energy is given by U, and m iS = C| s Ci iS . The elas- 
tic constant of the lattice is K, while the kinetic energy 
of lattice is neglected. The last term of Eq. (1) is the 
exchange coupling of two localized 1/2 spins (Sxi and 
St 2 ) to the spins of tt electrons at site il and i2, respec- 
tively. C\ is defined as (Cj^,Cj^) and a = {(T x ,a v ,a z ) 
are the Pauli matrices. We define a dimensionless cou- 
pling constant A = a 2 /{t$K) and dimensionless bond- 
lengths A, = a (<?,; - q i+ i)/t . 

We exactly diagonalizc the electronic and spin part of 
the Hamiltonian by the Lanczos algorithm. Therefore 
the obtained wavefunctions take account of all the corre- 
lation effects. The lattice deformation is treated classi- 
cally and optimized by means of the Hellmann-Feynman 
force equilibrium condition at zero temperature, 

A, A ^i.r ■ h.c.)Y (4) 

where (• • •} denotes the ground state expectation value. 
We perform these two procedures iteratively until the 
electronic state and lattice deformation converge. Since 
the model has the electron-hole symmetry, electron dop- 
ing gives the same results as hole doping. Only the results 
for hole doping will be discussed in the paper. 



First let us analyze the case of strong-correlation limit 
U/t — > oo. The half-filled Hubbard model is mapped 
onto the Heisenberg model with effective exchange inter- 
actions Jf s — —4t?/U (< 0), so that the electronic and 
spin parts of the present Hamiltonian is reduced to 

H eff = -E J ° S (SA+i-1/4) 

i 

- J(S« • S T i + S i2 • S T2 ). (5) 

The spin alignment predicted by the topological rule can 
be deduced from such a spin Hamiltonian in an intuitive 
way, if quantum fluctuations are neglected. Because of 
antiferromagnetic correlation between the neighboring 7r 
spins, the alignment between Sti and St2 can be classi- 
fied in terms of the parity of i2 — il: antiparallel if i2 — il 
is odd, but parallel if i2 — il is even, no matter whether 
the exchange coupling J is ferromagnetic or antiferro- 
magnetic. For J < 0, this intuitive picture is consistent 
with rigorous theorems—^i On the other hand, single 
hole doping induces Nagaoka ferromagnetic ordering in 
the limit of strong correlation^ Therefore the two local- 
ized spins favor parallel alignment. If U/t is finite, no 
simple rule exists for the doped w electron systems, and 
numerical methods are more effective way to study doped 
states. 

The model is investigated with the following param- 
eters in the unit of t (on the order of eV): U = 2.5, 
| J\ = 0.2 while A is allowed to vary from to 0.5. In the 
previous paper—, we assumed the exchange interaction 
much stronger (|J| = 1) than this value, and discussed 
mostly the case of ferromagnetic coupling. Here, weaker 
exchange interaction is applied to the model because the 
reduced value gives more realistic spin gap. Furthermore, 
J will be discussed for both ferromagnetic (F) J > and 
antiferromagnetic (AF) J < 0, although J has been sug- 
gested to be positive in the thianthrene-based molecular 
magnets^. 



III. RESULTS AND DISCUSSION 
A. Spin correlations and lattice deformations 

To investigate spin alignment with electron-lattice in- 
teraction, we first study a 10-site chain with A = 0.25 
(a typical value for polyacetylen o 33 ! 38 ) . The two local- 
ized spin are coupled to electron spin at both ends of the 
chain, i.e. il = 1 and i2 = 10. The spin correlations are 
plotted in Fig. El for both the half-filled (N e = 10) and 
the single hole doped (N e = 9) cases. 

The correlations between Sti and the other spins are 
shown in Fig. |2Ja) for the case of F coupling ( J > 0). The 
correlation between Sti and St2 is antiferromagnetic for 
the half-filled case, resulting in spin singlet 5 = 0. This 
antiparallel spin alignment can be understood based on 
the topological rule for the half-filled ground state. The 
spin correlation ((Sti • St 2 ) = —0.745) is near to —3/4 
of a singlet pair of the localized spins. 
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FIG. 2: Spin correlations (S T i-S;) for the half-filled (N e = N; 
dashed line) and the single hole doped cases (7V e = N— 1; solid 
line) with \J\ = 0.2, U = 2.5, A = 0.25, and N = 10. (a) F 
coupling (J > 0) and (b) AF coupling (J < 0). The inset in 
each figure shows schematically the change of spin alignment 
by hole doping. 



Positive correlation appears between the localized 
spins when the it electron moiety is doped by a hole. 
Similar spin alignment has been reported in the case of 
the stronger coupling— In the present case, this spin cor- 
relation ((Sti • St2) = 0.249) is near to 1/4 of a triplet 
pair of the localized spins. The spin correlation pattern 
in the right half of the chain (i > 6) is also reversed by 
the doping. Furthermore, the doped system turns out 
to be a spin quartet (5 = 3/2). The inset in Fig. [21(a) 
illustrates the alteration of spin alignment from the spin 
singlet to the spin quartet by such hole doping. This kind 
of low-spin to high-spin transition corresponds to the ob- 
servation in the aforementioned thianthrene derivative^ 

Figure H£b) shows the spin correlation in the case of 
AF coupling (J < 0). Antiferromagnetic correlation be- 
tween Sti and St2 (= —0.742) is found in the half- 
filled case with the total spin S = 0, which is also ex- 
pected from the topological rule. However, hole doping 
induces a spin doublet S = 1/2, while the correlation 
((Sti • St2) = 0.243) is ferromagnetic similar to the case 
of F coupling. The antiparallel alignment between the 
hole spin and the localized spins is responsible for the 
spin doublet as schematically shown in the inset of Fig. 

EDO- 

The lattice deformation patterns is plotted in Fig. |3| 
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FIG. 3: Lattice deformations A< for the half-filled (JV e = N; 
dashed line) and the single hole doped cases (N e = N — 1 ; 
solid line) with |J| = 0.2, U = 2.5, A = 0.25, and N = 10. (a) 
F coupling (J > 0) and (b) AF coupling (J < 0). 
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FIG. 4: (a) Spin correlations (Sti -St2) and (b) gap energies 
as functions of A for the half-filled (iV e = N; dashed line) 
and the doped (N e = N — 1; solid line) cases with J = 0.2, 
U = 2.5, A = 0.25, and N = 10. The inset in (b): the total 
spins of the ground and the first excited states. 



for both F and AF couplings. Quite similar profiles are 
obtained in both cases: half-filled case has bond alter- 
nation in the lattice, while the doped case contains very 
weak bond alternation especially around the middle of 
the lattice. The localization of the lattice relaxation is 
ascribed to a polaronic effect as well as a chain-end effect 
as shown in ref. 

We examine the effect of electron-lattice coupling on 
spin alignment by investigating spin correlation and gap 
energy between the ground state and the first excited 
state in the range of < A < 0.5. The A-dependences 
are plotted in Fig. [3] for F coupling. For both the half- 
filled and the doped cases, the spin correlation is barely 
affected by A. The inset in Fig. 01b) indicates that the 
lowest excitation takes place from singlet to triplet in the 
half-filled case, and from quartet to doublet in the doped 
case. In the half-filled case with A = 0.25, the spin gap 
(5.5 x 10 -4 to, about a few K) is comparable to the ef- 
fective couplings of thianthrene bis(nitronyl nitroxide) 7 . 
The spin gap gradually decreases as the electron-lattice 
coupling increases, because of the increasing tendency of 
local singlet-pair formation in it electrons. The gap in 
the doped case is an order of magnitude larger than that 
in the half-filled case, and decreases as A increases. Spin 
polarization is suppressed by strong bond alteration as 
electron-lattice coupling increases. The spin gaps of very 
weakly coupled systems tend to decrease with reduction 



of spin densities on the coupling sites. This A-dependencc 
of the spin gap is different from that reported in our pre- 
vious paper with a larger J—. In the case of strong ex- 
change interaction, the localized spins are firmly bound 
to the 7r-spins, and are hardly set free by low-energy exci- 
tations. Thus the spin gap would be mainly governed by 
excitations within the 7r-moiety. These 7r-excitation ener- 
gies are raised by the electron-lattice coupling, resulting 
in the increase of the spin gap with A in the stronger J 
case. 

For the AF coupling, the spin correlation and the gap 
energy are plotted as functions of A in Fig. |5] We see 
tendencies similar to the case of the F coupling, even 
though the manner of spin alignment in the doped state 
depends on the sign of the coupling, as already seen in 
the insets of Fig. [2J In the half-filled case with A = 
0.25, the spin gap (9.5 x 10 -4 to) is larger than that of 
the F coupling but is still fairly small against thermal 
fluctuations at room temperature. The suppression of 
electron-lattice coupling is desirable in molecular design 
with respect to the thermal stability. 



B. Effects of molecular structure 

Geometrical structures are very important to tt- 
conjugated molecular magnets. The influence on the 
doped state is different from that on the neutral one, as 
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pointed out in the previous papei***. It has been observec 
in the topological isomers of thianthrene bis(nitronyl ni 
troxide), that spin state depends on the positions of 
the radicals in the neutral case but has no position- 
dependence in the doped caset In this paper, we ex- 
amine two geometrical factors: the length of the chain N 
and the positions of the localized spins. 

First we examine the iV-dependence of spin alignment. 
Figure shows the correlations (Sti 4 St2) between the 
localized spins, the total spins, and the gap energies as 
functions of N in the case of the F coupling. To focus 
on the effect of the length, we maintain the localized 
spins at both ends of the chain. In the half-filled case, 
the topological rule governs the spin correlation and total 
spin perfectly. For even number of N, spin state is singlet 
(S = 0) with antiparallel alignment between Sti and 
St2 , while for the odd number case total spin is quartet 
(S — 3/2) with parallel alignment. 

Similar alternating behavior is also seen in the hole- 
doped systems, although the topological rule is generally 
only applicable for the ground state of the half-filled sys- 
tem. The hole doping changes the total spin and spin 
correlation (Sti ■ St2) as follows: 5 = 0^5 = 3/2 
with (S T i • S T 2) > — > (S T i • S T 2) < 0, if N is even; 
5 = 3/2^5 = with (S T i • S T2 ) < -> (S T i • S T2 ) > 0, 
if N is odd. It is noted that the spin correlation does not 
decay with increasing the chain length both in the doped 
and half-filled cases, while the spin gap decreases gradu- 



FIG. 6: (a) Spin correlations (Sti • St2), (b) total spins, and 
(c) gap energies as functions of N for the half-filled (N e = N; 
open circle) and the doped (N e = N — 1; filled circle) cases 
with J = 0.2, U = 2.5, and A = 0.25. 



ally as shown in Fig. EIc). In general, one-dimensional 
electron systems have no long-range charge (spin) order 
as a result of quantum fluctuations. The correlation be- 
tween a pair of electrons decays in a manner of power law 
with increasing distance. Decreasing electron correlation 
would weaken the interaction between the localized spins 
in our model, and reduce the spin gap as the chain be- 
comes longer. 

We turn to the case of the AF coupling shown in Fig. 
In the half-filled system, the alignment of Sti and St2 
is determined by the parity of N in the same manner of 
the F coupling: antiparallel if N is even, parallel if N 
is odd. In contrast to the F coupling, the total spin is 
singlet (5" = 0) for even N, and doublet (5" = 1/2) for 
odd N, as shown in Fig. Cfb). The single hole doping 
flips the localized spin from antiparallel to parallel with 
5 = 0^5=1/2 if N is even, and from parallel to 
antiparallel with S = 1/2 -> S = if N is odd. The 
./V-dependence of the spin correlation and of the spin gap 
is similar to the case of the F coupling. 

Next we proceed to the other geometrical factor, con- 
sidering the localized spins attached to inner sites of the 
chain. We pick up symmetric positions of the localized 
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spins as (il, i2) = (1, N), (2, N - 1), ■ • •, with TV = 10. 
In Fig. |S| the spin correlations (S T1 • S T2 ) , (S T1 • S c ) 

(S e = X)i=i ^i)) and the total spin are plotted as func- 
tions of il for the F coupling. In the half-filled case, these 
values are almost independent of il: the spin alignment 
is antiparallel with S = 0, as expected by the topological 
rule, since i2 — il is odd. The correlation between Sti 
and S c is negligible, because (S c ) is almost zero in the 
half-filled electron systems. 

For the doped case, the most interesting feature is that 
the relative alignment among the two localized spins and 
the hole spin depends on il and i2. In Fig. Ufa), the 
localized spins exhibit positive correlation (Sti • St2) ex- 
cept for the position (il,i2) = (4,7). Furthermore, these 
positively correlated spins can be classified into two cat- 
egories in terms of (Sti ■ S c ), namely the alignment be- 
tween the localized spin and the hole spin. They are 
parallel in the case of il = 1 , 3, and 5 but antiparallel 
in the case of il = 2 (Fig. Efb)). Correspondingly the 
doped systems are spin quartet (S = 3/2) for the for- 
mer three structures but spin doublet (S — 1/2) for the 
last one (Fig. Etc)). For (il,i2) = (4,7), the localized 
spins are in antiparallel alignment leading to spin doublet 
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(5 = 1/2). 

Figure|H|shows results for the AF coupling as the local- 
ized spins are attached to inner sites of the chain. Just 
like the F coupling, the half-filled systems do not ex- 
hibit the il-dependence: they are all spin singlet (S = 0) 
with antiparallel alignment between the localized spins. 
In the doped systems, the localized spins show parallel 
alignment in all the cases as shown in Fig. Ela). On 
the other hand, they show antiparallel alignment with 
the hole spin except for il = 2 (Fig. Efb)). This differ- 
ence in spin alignment results in the different total spin: 
Spin doublet (S = 1/2) for il = 1 ,3,4, and 5 and spin 
quartet (S = 3/2) for il = 2 (Fig. EJc)). 

The gap energies of the doped systems are presented 
in Fig. ^3 The gap is relatively large at il = 1 and 3 for 
both J > and J < 0, where the high-spin states appear 
in the case of J > 0. All the results above imply that 
molecular systems with particular structures are suitable 
for the purpose of spin-alignment control since they ex- 
hibit doping-induced high-spin and have large spin gaps. 



C. Dependence on J and U 

It is a key issue in the present study to understand 
the complex interplay of spin and charge degrees of free- 
dom in doped quantum systems. We have demonstrated 



C/3 

V 



A 



B9 

V 




M 

W 
& 

o 



60 

H 

<D 

W 
& 

Co 

o 




FIG. 9: (a) Spin correlations (Sti-St2) and (b) (STi-S e ), and 
(c) total spins as functions of il, with i2 = N -\- 1 — il, for the 
half-filled (JV e = N; open circle) and the doped (N e = N - 1; 
filled circle) cases with J = —0.2, U = 2.5, A = 0.25, and 
N = 10. 



the spin-alignment control by hole-doping using fixed pa- 
rameters of the on-site Coulomb repulsion U and the 
exchange interaction J. However, a subtle balance be- 
tween these interactions governs the alignment between 
the separated radical spins. It leads to a question: Does 
spin alignment change dramatically with varying the 
strengths of U and J, or just depend on them quanti- 
tatively? As already shown in Fig. the spin gap is 
sometimes very small in the doped systems, hence the 
order of the present lowest excited and ground states 
could be interchanged for another set of parameters. Fur- 
thermore, the spin alignment depends on the topological 
structure as shown in Sec. IIII Bl Thus, it is necessary to 
investigate the dependence on U and J for various molec- 
ular structures. In the following discussion, the electron- 
lattice coupling is omitted as it plays less important role 
than the electron-electron interaction. 

In the case of F coupling, the phase diagrams in the 
parameter space of J and U are shown for il — 2, 4, and 
5 in Fig. The doped systems with il — 1 or 3 are al- 
ways spin quartet in the calculated range of < J < 0.5 
and < U < 8. As the on-site Coulomb interaction in- 
creases over a critical value about U ~ 0.4 for il = 2, or 
U ~ 2.0 for il — 4, the systems are transferred from the 
spin quartet to the spin doublet, and the phase bound- 
aries are almost independent of the exchange interaction. 



FIG. 10: Gap energies of the doped system: (a) J > and 
(b) J < as functions of il (i2 = N+ 1 - il), with \ J\ = 0.2, 
U = 2.5, A = 0.25, and TV = 10. 



The system with il = 5 has much strong dependence on J 
and U. It also undergoes a transition from the spin quar- 
tet to the spin doublet, but occurs in stronger Coulomb 
interaction U > 4.0. Detailed analysis reveals that the 
spin gap goes continuously to zero when the system ap- 
proaches to the phase boundary. 

Similar transitions are also found in the case of AF 
coupling. Figure 1121 shows the phase diagrams for the 
systems with il = 2 and 4, which are transferred from 
spin doublet to quartet at U ~ 0.4 for il = 2 or at 
U ~ 2.0 for il = 4 in a way opposite to the case of F 
coupling. The transitions occur at similar values of U in 
both the F and AF couplings. However, J-dependence is 
opposite between them. The systems with il = 1,3, and 
5 do not exhibit phase transitions in the calculated range 
of -0.5 < J < and < U < 8. 

To understand the origin of these phase transitions, 
we calculated the doped states in the Id Hubbard model 
(N = 10, iVf = 5, Ni = 4), which is nothing but setting 
J = in our model. Their spin densities are plotted 
as functions of site(i) in Fig. I13f a). The spin density is 
positive at every site if the on-site Coulomb interaction is 
absent or quite weak. As the strength of U increases, spin 
densities become negative first at i — 2 and 9 for U > 0.4, 
and then at site i = 4 and 7 for U > 2.0. The spin 
densities at the other sites remain positive irrespective 
of the strength of U. We also plot the spin correlation 
(Sn ■ S i2 ) with i2 = N + 1 - il in Fig. [H». The 
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FIG. 11: Phase diagrams of the doped system (7V e = TV — 1) 
for (a) il = 2, (b) il = 4, and (c) il = 5, (i2 = N + 1 - il) 
with J > 0, A = 0, and N = 10. 



correlations are positive for il = 1 and 3, negative for 
il = 4 and 5, and almost zero for il = 2. 

Returning to our original model with the local- 
ized spins attached, we can understand the position- 
dependent alignment of the spins shown in Figs. |3]and[!5] 
from the spin densities of the it electrons shown in Fig. 
ITST a). In the case of F coupling, each localized spin tends 
to be in the parallel alignment with the 7r-spin at its cou- 
pling site. If the coupling site has a positive spin density, 
the localized spin is in the parallel alignment with the 
hole spin leading to a spin quartet, as illustrated in Fig. 
I14f a'l. On the other hand, the radical spin is in the an- 
tiparallel alignment with the hole spin if the coupling site 
has a negative spin density, resulting in a spin doublet as 
schematically shown in Fig. I14f b1. In addition, we can 
also explain the critical values of U in the phase diagrams 
of Fig. ^2 in the limit of J — > 0: The transition point is 
close to the value of U where the spin density becomes 
negative in Fig. I13f a~). The spin alignment in the case of 
AF coupling can also be understood in a similar way as 
shown in Fig. I14f c) and (d). This picture is consistent 
with the spin correlation (S,i • S^}, which is positive for 
large U at il = 1, 2, and 3 as shown in Fig. I13f b). As for 
il = 5, the spin alignment for both F and AF couplings 
can be explained based on the spin densities, although 
the spin correlation (S;i ■ S;2) is negative. 

The absolute value of the spin density at il = 4 is 
very small, and the parameter set is close to the phase 
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FIG. 12: Phase diagrams of doped system (N e = N — 1) for 
(a)il = 2 and (b)il = 4,(i2 = N + 1 - il) with J < 0, A = 0.0 
and N = 10. 



boundary as shown in Figs. I12f b). The sign of spin den- 
sity at this site can be changed to positive by a finite 
electron-lattice coupling. As a result, spin doublet ap- 
pears in the case of AF coupling, as schematically shown 
in Fig. I14f c). The spin correlation (S4 • S7) is negative 
as shown in Fig. 1131^ ). and takes control of the spin 
alignment in the case of F coupling. This spin alignment 
leads to antiferromagnetic correlation between the local- 
ized spins, and results in a spin doublet. Such situation 
appears also at il = 5 for the case of stronger coupling 
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in our previous paper 



IV. SUMMARY 

In summary, the control of spin alignment by charge 
doping is studied in the theoretical model of polyene- 
based molecular magnets. By mean of the exact diag- 
onalization technique, we study the spin alignment, the 
total spins and gap energies in the doped and the half- 
filled states. 

Low-spin to high-spin transition is realized in such 
polyene-based molecular magnets by one-hole doping 
into the 7r-conjugated moiety. The doping-induced spin 
alignment depends on the molecular structure and the 
strengths of the interactions. Alternating behavior of 
spin alignment is demonstrated in the hole-doped sys- 
tems with growth of the chain length. Variation of spin- 
coupling positions leads to different alignment among the 
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spin alignment can be understood in terms of spin den- 
sity and correlation of n electrons, in contrast to the it 
topological rule for the half-filled case where antiferro- 
magnetic correlation in tt electrons governs spin align- 
ment. 

We have shown that charge doping is an effective way 
to realize controllable spin alignment in molecular mag- 
nets. As an important steD towards molecular designing. 
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FIG. 13: (a) Spin densities (n^f — 7ii^)/2 as a function of i; (b 
spin correlations (Sji -Sia) as a function of il (i2 = N+l — il 
for the doped state in the f d Hubbard model (N = f 0, Nf = 
5, iVj = 4) with open boundary condition. The dash-do1 
dash, and solid lines correspond to U = O.f, 0.4, and 2.C 
respectively. 



localized spins and the hole spin. To understand behav 
iors in doped system, we investigate the dependence c 
spin alignment on J and U in detail. The antiferromag 
netic and ferromagnetic exchange interactions have dii 
ferent effect on doped molecules. The overall behavior c 
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FIG. f4: Schematic picture of intramolecular spin alignment 
of doped system. In the case of J > 0: (a) S — 3/2 and (b) 
S = 1/2; in the case of J < 0: (c) S = 1/2 and (d) S = 3/2. 
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